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Abstract 

Abstract: We provide a Hamiltonian analysis of the Mixmaster 
Universe dynamics showing the covariant nature of its chaotic behavior 
with respect to any choice of time variable. We construct the appro- 
priate invariant measure for the system (which relies on the existence 
of an "energy-like" constant of motion) without fixing the time gauge, 
i.e. the corresponding lapse function. The key point in our analysis 
consists of introducing generic Misner-Chitre-like variables containing 
an arbitrary function, whose specification allows one to set up the 
same dynamical scheme in any time gauge. 



PACS number (s): 04.20. Jb, 98.80.Dr 



The original Belinski-Khalatnikov-Lifshitz (BKL) analysis of the asymp- 
totic dynamics of Bianchi type VIII and IX cosmological models (the 
so-called Mixmaster universes 0) showed the existence of a chaotic behavior 
in their approach to the initial cosmological singularity. 

Probably the most suitable description of this BKL oscillatory regime, 
especially in view of its chaotic nature, is via the Hamiltonian formulation of 
the Mixmaster Universe dynamics expressed in Misner-Chitre-like variables 
The advantage of this approach, apart from the immediate inter- 
pretation of the system evolution as a chaotic scattering process, consists of 
replacing the discrete BKL map by a geodesic flow in a space of continuous 
variables. In this approach a key result is obtaining an invariant measure for 
the Mixmaster evolution (see also [HIIZI)- 

In spite of the achievements made in describing the Mixmaster stochas- 
ticity, however, because of its relativistic nature, the dynamics can be viewed 
from different reference frames, thus leaving open the question of the covari- 
ance of the observed chaos. 

Interest in these covariance aspects has increased in recent years in view 
of the contradictory and often dubious results that have emerged on this 
topic. The confusion which arises regarding the effect of a change of the 
time variable in this problem depends on some special properties of the Mix- 
master model when represented as a dynamical system, in particular the 
vanishing of the Hamiltonian and its non-positive definite kinetic terms (a 
typical feature of a gravitational system). These special features prevent the 
direct application of the most common criteria provided by the theory of 
dynamical systems for characterizing chaotic behavior (for a review, see |H]). 

Although a whole line of research opened up, following this problem of 
covariant characterization for the chaos in the Mixmaster model Pl^CSl; the 
first widely accepted indications in favor of covariance were derived by a frac- 
tal formalism in ^3] (see also jlSI)- Indeed the requirement of a covariant 
description of the Mixmaster chaoticity when viewed in terms of continu- 
ous dynamical variables, due to the discrete nature of the fractal approach, 
leaves this subtle question open and prevents a general consensus from being 
reached. 

Here we show how the derivation of an invariant measure for the Mix- 
master model (performed in IZ| within the framework of the statistical 
mechanics) can be extended to a generic time gauge (more directly than in 
previous approaches relying on fractal methods jEj) provided that suitable 
Misner-Chitre-like variables are chosen. More precisely, after a standard 
Arnowitt-Deser-Misner (ADM) reduction of the variational problem when 
written in terms of a generic time variable, we show how asymptotically close 
to the cosmological singularity, the Mixmaster dynamics can be modeled by 
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a two-dimensional point-universe randomizing in a closed (in the sense below 
precised) domain with fixed "energy" (just the ADM kinetic energy); since 
it is natural to represent such a system by a microcanonical ensemble, then 
the corresponding invariant measure is induced by the Liouville one. The 
key point, making our result reliable, is that we consider an approximation 
dynamically induced by the asymptotic vanishing of the metric determinant. 

As well-known ^H]; the Mixmaster dynamics is governed by the La- 
grangian (written in the Misner variables j2]) ^H] 
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where (-I-) and (— ) refer respectively to Bianchi type VIII and IX, and the 
anisotropy parameters Hi (i = 1, 2, 3) are the functions j3] 
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In the limit D ^ corresponding to approaching an initial singularity, the 
second set of three terms of the above potential are clearly negligible with 
respect to the first set of three terms (excluding the particular cases when 
two or three anisotropy parameters Hi coincide). 

Next introduce generic Misner- Chitre-like variables by the transformation 
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^2 - 1 sin 9, 



where <oo,0<^^<27r and r(r) stands for a generic function of r 

(taken to be identical to r by Chitre, who also sets ^ = cosh^, < ( < 00, 
see jni)- In terms of these variables the Lagrangian becomes 
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where D is 

D = exp{-3^e^(")}. (6) 

Hence, by the usual Legendre transformation we obtain the conjugate mo- 
menta 

12D / p drV , 

and reformulating the dynamical problem in terms of the corresponding 
Hamiltonian variational principle leads to 
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S J [p^e + Pee' + PrT - -^^nj dr^ = 0, (8) 
where 

n = --g, + p,^ {e-l)+^^+ 24Fe-. (9) 

By varying (jH)) with respect to we get the constraint H = 0, whose 
solution is 

-p. ^ ^Hadm = ^Ve^ + 24Ve^^, (10) 
dr dr 



where 
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Thus we reduce the variational principle (jH)) to the simpler one 



S J (Pc^' + Pe0' - T'Hadm) dr] = 0; (12) 



associated with this elimination, and from the first of ((Zj) and |T(H) follows 
the time gauge relation 

, , 12D yrdT , , , 

^{i) = W, e^^-r'. (13) 

riADM dT 

Although now we choose the natural time gauge r' = 1 and rewrite the 
variational principle in terms of the time variable F (indeed the Hamiltonian 
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equations are equivalently viewed through the two time variables F and r) 
leading to the variational principle 



nevertheless for any choice of time variable r (i.e. t = rj), there exists a 
corresponding function F (r) (i.e. a set of Misner-Chitre-like variables leading 
to the scheme (^H)) defined by the (invertible) relation 

As a consequence of the variational principle (fT^ we easily get the funda- 
mental relation 

When expressed in terms of the Misner-Chitre-like variables, the anisotropy 
parameters Hi [i = 1,2,3) are independent of the variable F 



Hi = ^- ^ (cos g + vising 



H2 = ^- "^^g'^ cosg- vising) (17) 
H3 = - + 2 ^^"~^ -cosg, 

so as D vanishes asymptotically approaching the initial singularity, the po- 
tential term can be modeled by the time- independent (see FiglH). potential 
walls 

u^ = eoo{Hi{^,e)) + e^{H2{^,e)) + e^{Hs{^,e)) , (is) 

©oo (x) = +00 if X < 0, 
©oo (x) = if X > 0, 

and therefore in the region FIj^ where the potential vanishes, we have by (fT?)|l 
dH-ADM /dT = 0, i.e. e = E = const (by other words the ADM Hamiltonian 
asymptotically approaches an integral of the motion). 

The invariant measure for this system can then be derived following the 
same approach presented in [3 [7j . 

Indeed the Mixmaster dynamics in this approximation is equivalent to 
a two-dimensional point- universe moving within closed potential walls in a 
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surface of constant negative curvature (the Lobachevsky plane [5|), described 
by the hne element 



where is a constant. From a geometrical point of view, the domain defined 
by the potential walls is not strictly closed, since there are three directions 
corresponding to the three corners in the traditional ADM-Misner picture 
from which the point- universe could in principle escape (see Fig^. 

However, as well known, for the Bianchi models under consideration, the 
only case in which there exists an asymptotic solution of the field equations 
with this behavior corresponds to two scale factors being equal to each other; 
but as shown in jT^j, these cases are dynamically unstable and therefore 
correspond to sets of zero measure in the space of the initial conditions. For 
this reason it does not make sense to speak of a probability to reach certain 
configurations and the domain is de facto dynamically closed. 

The bounces (billiard configuration) against the potential walls together 
with the instability of the geodesic flow on a closed domain of the Lobachevsky 
plane cause the point-universe dynamics to have a stochastic feature. In- 
deed just the condition that the potential-free motion takes place in the 
Lobachevsky plane is itself not sufficient to infer the chaoticity of the dy- 
namics, which is true for any Bianchi type model. It is the compactness of 
the domain bounded by the potential walls that guarantees that the geodesic 
instability is transformed into a real stochastic behavior. On the other hand, 
the possibility to speak of a stochastic scattering is justified by the constant 
negative curvature of the Lobachevsky plane and therefore these two notions 
(compactness and curvature) are both necessary for our considerations. 

From the statistical mechanics point of view, such a system admits an 
"energy-like" constant of motion which corresponds to the kinetic part of 
the ADM Hamiltonian e = E. Since the point-universe randomizes within 
the phase-space associated with the closed domain Uh, its stochastic dynam- 
ics admits a suitable ensemble representation which, in view of the existence 
of the "energy-like" constant of motion, has to have the natural feature of a 
microcanonical one. Therefore the stochasticity of this system is appropri- 
ately described in terms of the Liouville invariant measure 



where S (x) denotes the Dirac function, characterizing the microcanonical 
ensemble. 

It is important to note that the particular value taken by the constant e 
{e = E) cannot influence the stochastic property of the system and must be 
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fixed by the initial conditions. To remove this useless information from the 
statistical dynamics under consideration, we must integrate over all admissi- 
ble values of e. To do this it is convenient to introduce the natural variables 
(e, 0) in place of {p^,pe) by 

Pf = , ,f cos 6 



Pe = e^e-lsin^, (21) 

where < < 27r. Integrating over all possible values of e removes the Dirac 
delta function, leading to the uniform (normalized) invariant measure 

dfi = didOdcp^. (22) 

It is important to emphasize that the approximation on which our analysis 
is based (i.e. the potential wall model) is reliable since it is dynamically 
induced no matter what time variable r is used. 

The analysis developed in |7J shows that it is possible to derive asymp- 
totically the Liouville theorem for the reduced phase space {^,6',0}. The 
relative continuity equation for the probability distribution w (F, 9, 0) can be 
expressed as 

dw d^ dw dO dw dcf) dw 

sin0 dw ^smcj) dw 



One can immediately check that the invariant measure ()22|) is a stationary 
solution of this continuity equation. 

In fact, by applying the Landau- Raichoudhury theoren:0 near the initial 
singularity (which occurs by convention at T = 0, where T denotes the 
synchronous time, i.e. dT = —N (r) dr), we easily obtain the result that D 
vanishes monotonically (i.e. for T — we have d In D/dT > 0). In terms of 
the adopted time variable r (D — > ^ F (r) ^ oo), we have 

dlnD dlnDdT dlnD , , 



^This theorem, based on the mathematical assumptions underlying the dynamics of 
the Einstein equations, states that in a synchronous reference frame there always exists a 
value of the time at which the metric determinant vanishes and that in this time variable 
the zero is approached monotonically jl8| . 
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and therefore D vanishes monotonically for increasing r as soon as rfF/dr > 
according to (fT3|) . 

The key point of our analysis is that any stationary solution of the Liou- 
ville theorenu, like ()2()|1 . remains valid for any choice of the time variable r; 
indeed the construction presented in ^ of the (restricted) Liouville theorem 
shows the existence of such properties even for the invariant measure 
, Clearly the knowledge of the invariant measure (j^^ provides a satisfac- 
tory statistical representation of the system for any choice of time variable, 
since it allows one to calculate the asymptotic average values (as well as 
higher order moments) of any dynamical variable involved in the problem. 

We conclude by remarking that when approaching the singularity F (r) — ^ 
oo (i.e. TCadm E), the time gauge relation (fT3j) simplifies to 

^.Mir)^^e-^r, (25) 

According to the analysis presented in 0, the asymptotic functions ^ (F) , 9 (F) , (F) 
during free geodesic motion are governed by the equations 
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Once the solution C, (F) is obtained in the parametric form 
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r(0) 

p = \J a? + sin^ (j) a,b = const. E ^ (27) 

equation (f^3j) for free geodesic motion reduces to a simple differential equa- 
tion for the function F (r) . 

However, as shown by our analysis, the global behavior of C, along the 
whole geodesic flow is described by the invariant measure p2|l and therefore 
relation indeed takes a stochastic character. In other words, if we assign 

^When using a generic time variable r, the right-hand side of equation 1231 in general is 
no longer (asymptotically) vanishing, but yet negligible with respect to the left-hand one. 
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one of the two functions F (r) or (r) to an arbitrary analytic functional 
form, then the other one will exhibit stochastic behavior. Finally, by retain- 
ing the same dynamical scheme adopted to construct the invariant measure, 
we see how the one-to-one correspondence between any lapse function (r) 
and the associated set of Misner-Chite-like variables guarantees the co- 
variant nature with respect to the time-gauge of the Mixmaster universe 
stochastic behavior, when viewed in the framework of statistical mechanics. 

Of course, our analysis leaves open the question concerning the covariance 
of the Mixmaster chaos with respect to the choice of configurational variables 
which are not Misner-Chitre-like. 



Acknowledgments 

We are very grateful to Vladimir Belinski and Remo Ruffini for their valuable 
comments on this subject. Robert Jantzen and Carlo Bianco are thanked 
respectively for the help in rewriting this manuscript and in producing the 
figure. 

References 

[1] Belinski V A, Khalatnikov I M and Lifshitz E M (1970) Adv. Phys., 19, 525. 

[2] Misner C W (1969) Phys. Rev. Lett, 22, 1071. 

[3] Chitre D M (1972) Ph.D. Thesis, University of Maryland. 

[4] Chernoff D F and Barrow J D (1983) Phys. Rev. Lett, 50, 134. 

[5] Kirillov A A and Montani G (1997) Phys. Rev. D, 56, n. 10, 6225. 

[6] Barrow J D (1982) Phys. Rep., 85, 1. 

[7] Montani G (2000) in: The Chaotic Universe, Gurzadyan V G and Ruffini R 
(eds), 173, World Scientific, Singapore. 

[8] Hobill D, Burd A and Coley A, (eds) (1994) Deterministic Chaos in General 
Relativity, World Scientific, Singapore. 

[9] Pullin J, (1991) Relativity and Gravitation: Classical and Quantum, Proc. 
Vnth SILARG Symp. Cocoyoc Mexico 1990, J C D'Olivio et al. (eds) 

[10] Rugh S E in (1994) Deterministic Chaos in General Relativity, Hobill D et 
al. (eds), World Scientific, Singapore. 



8 



[11] Biesiada M (1995) Class. Quantum Grav. 12, 715. 

[12] Szydlowski M (1993) Phys. Lett. A , 176, 22. 

[13] Szydlowski M and Szczesny J, (1994) Phys. Rev. D 50, 819. 

[14] Cornish N J and Levin J J (1997) Phys. Rev. Lett, 78, 998; (1997) Phys. 
Rev. D, 55, 7489. 

[15] Motter A E and Letelier P S, preprint |gr-qc/ 001100T 

[16] Misner C W, Thorne K S and Wheeler J A, Gravitation, ed.by W H Freeman 
and C, New York (1973) p. 812. 

[17] BeUnski V A and Khalatnikov I M (1969) Sov. Phys. JETP, 29, 911. 

[18] Landau L D and Lifshitz E M, "Classical Theory of Fields", (Mir), (1975), 
4th Edition, Addison- Wesley, New York. 



9 



^ 10 



Figure 1: Asymptotic potential domain Fh in the Lobachevsky plane 9,^. 
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